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SUMMARY 


This document describes the computational aspects of propeller noise prediction 
in the time domain and the use of the high speed propeller noise prediction pro- 
gram ASSPIN (Advanced Subsonic and Supersonic Propeller Induced Noise). The 
code, which was formerly called the DFP-ATP (Dunn-Farassat-Padula Advanced 
Technology Propeller) noise prediction program, was developed at NASA Langley Re- 
search Center (LaRC) by the authors and F. Farassat and S. Padula of NASA LaRC. 
ASSPIN is the latest version of a sequence of propeller noise prediction codes that are 
based on the theoretical formulations of Farassat. These formulations are valid in both 
the near and far fields. Two such formulations are utilized by ASSPIN. One formulation 
is used for subsonic portions of the propeller blade, and the other for the transonic and 
supersonic regions. Switching between the formulations is performed automatically. 
ASSPIN incorporates advanced blade geometry and surface pressure modelling, adap- 
tive observer time grid strategies, and relative to the previous version of the propeller 
noise prediction code, contains enhanced numerical algorithms that result in reduced 
computational time. In addition, the ability to treat the nonaxial inflow case has been 
included. 

INTRODUCTION 

Advanced Technology Propellers (ATP) or Propfans are a new generation of 
propellers that feature highly loaded blades that are thin and swept back (Figure 1). 
Propfan propelled aircraft are able to attain the same high speed, high altitude flight 
as current jet technology, but with a 20 to 25 percent reduction in fuel consumption 
(reference 1). Counter-rotating propfans (Figure 2) reduce swirl losses in the slipstream 
providing additional reduction in fuel consumption. 

Producing blade designs that suppress propeller noise is an active area of research 
in the field of aeroacoustics. The construction and acoustical testing of new propeller 
designs is expensive and time consuming. Analytical methods, on the other hand, 
provide the engineer with a cost effective alternative to the experimental determination 
of propeller noise. This work is motivated by the need for an accurate analytical 
propeller noise prediction method that includes all computationally feasible physical 
effects. 

ASSPIN is a computer program that predicts the noise generated by propellers 
operating at subsonic, transonic, or supersonic helical tip speeds in either single-rotation 
or counter-rotation mode. The prediction method is based on two theoretical time 



Figure 1. Single- Rotation Propeller Figure 2. Counter-Rotating Propellers 

domain formulations of Farassat (references 2,3). One formulation is valid for subsonic 
regions of a propeller blade, and the other, though valid for all speeds, is used only for 
transonic and supersonic portions of the blade. The formulations, without modification, 
apply in both the near and far fields, and to the nonaxial inflow case. Input to ASSPIN 
consists of blade kinematic and geometric quantities and either steady or periodically 
unsteady blade surface pressure. Program output includes propeller power and thrust 
coefficients, and the periodic acoustic signatures and acoustic spectra for the thickness, 
loading, and combined noise. 

ASSPIN is the latest version of a sequence of noise prediction codes (references 4,5,6) 
that are based on the theoretical formulations of Farassat. Relative to the previous 
version of the code (reference 6), ASSPIN contains enhanced numerical and analytical 
methods that are designed to reduce computational time and increase accuracy. These 
facets of ASSPIN are explored in the Computational Strategy section. 

ASSPIN has been used as part of a comprehensive propeller analysis package 
that includes state-of-the-art methods in propeller aeroacoustics, aerodynamics, and 
aeroelasticity. It has been demonstrated that acoustic results generated by ASSPIN 
within the propeller analysis code agree well with actual flight data (reference 7). 


LIST OF SYMBOLS AND ABBREVIATIONS 

Aj, A 2 reference frame rotation matrices 


\Mt + Ai <i ; b = 


B l 


components of b in direction of blade principal curvatures, i — 1,2 


2 


bn 

b * n 

CRP 

counter-rotation propeller 

C Q 

speed of sound in undisturbed fluid 

c P 

propeller power coefficient 

C T 

propeller thrust coefficient 

dB 

decibels 

dS 

element of blade surface area 

5i 

li 

o 

equation describing blade surface in reference frame fixed to undis 
turbed medium 

F{y\ x, t ) 

o 

0 

1 

*+o 

H 

local mean curvature of blade surface 

h n 

A M n + A] cos 6 

k(y, t ) 

function defined so that the intersection of the surfaces k = 0 and 
/ — 0 specifies the bounding curves of an open piece of the blade 
surface; intersection of the region k > 0 and the surface / = 0 
defines an open piece of the blade surface 

K (y; x, t) 

o 

0 

1 

-to 

M 

local Mach vector based on c 0 ; M = M 

Mau 

M • z7, M with respect to a medium fixed observer 

M, Mi 

time derivative of M 

Mn 

M • n 

M r 

M • r 

Mp 

jDrojection of M onto local normal plane to blade surface edges; 

M p = M p 

— * 


M t 

projection of M onto local tangent plane to blade surface; Mt = 

M t 

n, n-i 

unit normal to blade surface 

n b 

number of blades 


3 



p 

p 

Pb 

PC A 
Q 

r, ri 
r, ri 

r P 


SPL 

SRP 

t 

h 

v F 

V 

x 


y 

Y 


acoustic pressure 

blade surface pressure expressed in medium fixed coordinates 
time derivative of p 

blade surface pressure expressed in blade fixed coordinates 
pitch change axis 

relative chord curvilinear coordinate 
x -y,r = ( t^I 
r/r 
r • V 

unit vector in the direction of the projection of r onto the local 
normal plane to the surface edges 

sound pressure level 

single-rotation propeller 

observer time 

projection of r onto local tangent plane to blade surface 
forward flight velocity 
velocity vector 

observer position expressed in medium fixed coordinates 
source position expressed in medium fixed coordinates 
source position expressed in aircraft fixed coordinates 


Greek Symbols 

o inflow angle 

d'y elemental arc length along acoustic planform edges 

dr elemental arc length on curve formed by the intersection of the 

collapsing sphere and the moving blade for a fixed source time 

8 Dirac delta function 

d£ element of acoustic planform area 
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e 

v 

m 

A 

A 

A 0 

P 

v 


P 0 
°b 

< 711 , <722 

T 


angle between n and r 

blade surface point expressed in blade fixed coordinates 
spanwise curvilinear coordinate 
yj 1 + Ml - 2M n cos Q 
a/a 2 + sin 2 0 

yj M 2 cos 2 0 + ^1 — Mp • r p sin 

components of M* in direction of principal curvatures, i = 1,2 

local inward unit geodesic normal to the curve given by the inter- 
section of the surfaces k = 0 and / — 0 

density of undisturbed fluid 

length parameter along b on / — 0 

— * — » — — * 

two components of the tensor t{t\ — Mt\ 
source time 


ip 

n, 


K l, «2 


K m, K t, K b 


A 

Ai 

X 


angle between r and edge of surface given by the intersection of the 
surfaces k = 0 and / = 0 

n x cJ 

local principal curvatures of the surface / = 0 

local normal curvatures of the surface / = 0 in the direction of Mt, 
< 1 , and 6, respectively 
cos 9 — M n 

A 2 

cos 6 + M n 
A 2 

observer position expressed in blade fixed coordinates 


uj angular velocity; u> = \lj\ 

THEORETICAL FORMULATIONS 


Small perturbations in fluid quantities due to the motion of a body through an 
otherwise quiescent, ideal fluid are governed by the equations of linearized acoustics. 
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Therefore, in the region exterior to the body, the acoustic pressure, p'(x, t), is a solution 
of the homogeneous wave equation subject to the boundary condition that the fluid 
cannot penetrate the surface of the booty. 

Using generalized function theory, the restricted domain of the above problem, be., 
the body’s exterior, can be extended to include the entirety of space by the process 
of embedding (reference 8). Applying this procedure to the equations of linearized 
acoustics leads to the Ffowcs Williams- Hawkings equation without the quadrupole term 
(reference 9): 

1 ~d^ * 

where 

and the bars over the differential operators denote generalized differentiation (reference 
8). The boundary condition and the pressure on the surface of the body appear as 
source terms in (1). The advantage of the embedding technique is that the Green’s 
function for the wave equation in unbounded space is known, thus the solution of (1) 
can be obtained analytically provided the aerodynamic pressure on the surface of the 
body is known. 

There are several different ways to express the solution of (1). Two forms of the 
solution are particularly effective for numerical evaluation and are considered here. The 
first form, known as Formulation 1A (reference 2), is valid over source regions that 
move toward the observer point, x, at subsonic speed. The second form of the solution 
is referred to as Formulation 3 (reference 3). Formulation 3 is valid for the entire 
body regardless of its speed, but for computational reasons that are discussed later in 
this section it is used only for supersonic and transonic portions of the body. In the 
subsequent discussions, the acoustic pressure will be written as 

p(x, t) = p T (x, t)+p' L (x, t ), 

where p' T and p' L are referred to as the thickness and loading noise components, 
respectively. 

For Formulation 1A, the acoustic pressure components are written in terms of 
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surface integrals over the moving body as follows: 


4tt p' (x, t) = c 0 p 0 


J 

/= 0 


Mn (rMifi + c 0 M r c 0 M‘ ) 


r 2 (1 - M r ) 


dS, 


(2a) 


J ret 


and 


47 rp 




/=0 


p cos 8 

r [p(cos«-M„)| 

r ( 1 — Mr ) 2 

f l 0 [ rHl - M r y \ 


dS 


ref 


+ 


i / 

C 0 J 
/=0 


p cos # | rMifi + c 0 M r — c 0 M 2 ^ 


" 2 (1 - M r )' 


dS. 


(2b) 


J rei 


The subscript T€ t appearing in the above formulae implies that the integrands are to be 
evaluated at the retarded (or emission) time t — r/c 0 . Observe that evaluation of (2a) 
and (2b) involves the calculation of kinematic, geometric, and aerodynamic quantities 
on the surface of the moving body. In particular, note that the geomctiic quantities 
required to evaluate (2a) and (2b) depend on the mathematical surface representation 
(viz. r, M r , r,*, and / = 0) and it’s first derivatives (viz. M n , 0, and dS). 


Portions of the body that move toward the observer at tiansonic or supersonic 
speed require a different formulation than the above. This is because the Doppler 
factor, 1 - M r , appearing in (2a) and (2b) becomes small at transonic speeds relative 
to the observer, rendering the integrals improper. Consequently, the analysis upon 
which these equations are based is invalid. With minor modifications, (2a) and (2b) 
are valid for strictly supersonic regions of the moving body. However, experience has 
shown that it is computationally difficult to differentiate between nearly transonic 
and supersonic source regions. These considerations prompted the development of 
Formulation 3 which, as previously mentioned, is mathematically valid for all flight 
regimes. In contrast to Formulation 1A, Formulation 3 is written in terms of surface 
integrals over open portions of the acoustic planform surface and line integrals over the 
edges of the open portions. The acoustic planform is generated by open portions of the 
moving surface (a panel, for example). Using Formulation 3, the thickness and loading 
components of the acoustic pressure can be written as 
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4vr I /- 


p o c o 


p' T (x, t) = J i 


F=0 

AT>0 




dE + 


ret F=0 

K >0 


/; 


F=0 
K = 0 


A/2g„ + M n M< 


at/ 


1 M%Q f + Q' f + Q 
A 

d-l, 


dZ 


ret 


(3a) 


ret 


and 


4?r p' L (x,t)= J ~ 


F = 0 
K>0 


pQ' 


N 


A 


re/ 


r i [K?f + 7-Pb - 

dZ+ j II da b 


A 


F=0 

/O0 




J re/ 


f I 

pQe 

J r 

’=0 

A o . 




re/ 


A'zrO 


respectively, where 


Q' n = A [2Ai (cos 0 - M n ) + 1] , 


1 


Qf~ (-A ~ 9 ) Mn + 2b 2 + Kian + «2 <7 22 ~ 2Hh r 

Co 


+ —n ■ 

c„ 


A2 

Mt ~ h 


1 


<4 = 2Mn 


A 2 - 2A6+ tX? + 2AA llr 


- (a Mn -n-b)+ Kut l B l + K 2 fl 2 B 2 

C n \ / 


Q" f = — (4 - ft • M t ) + K m M 2 - 2 HM 2 , 

C n ' ' 


and 


(3b) 


Q e — AM ai/ A - AjTj/. 


In equations (3a) and (3b), the notation F — 0, A >0, defines an open piece of the 
acoustic planform surface, and F = 0, A =0 defines its boundary. The quantities F 
and K are defined by 


F(y; x,t) = f (y, t - r/c 0 ) and K ( y ; x, t ) = k (y, t - r/c 0 ) , 
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with / = 0, k > 0 defining an open portion of the moving surface and / = 0, k = 0 its 
boundary. 

The acoustic planform (or E surface) is a surface in space generated by the 
intersection of the surface of the moving body and a sphere, centered at x , whose 
radius collapses at the speed of sound to the point x at the time t. In practice, the 
surface integrals of (3a) and (3b) are best calculated by the so called collapsing sphere 
method (reference 10), in which the relation 

dE c 0 cWdT /.x 

A sin# 

is utilized. In (4), the source time integration is over the period(s) of time in which the 
collapsing sphere, given by the equation 


y| 


= t — r 


( 5 ) 


and the moving body, represented by /(y, r) = 0, intersect. For fixed r, dr denotes 
elemental arc length along the curve) s) of intersection of the collapsing sphere and the 
moving body. 

Similar to the integrands of Formulation 1 A, the evaluation of (3a) and (3b) involves 
the calculation of geometric, kinematic, and aerodynamic quantities. However, the 
integrands of Formulation 3 contain many more terms than those of (2a) and (2b). 
Furthermore, second derivatives of the mathematical surface representation must be 
computed. This is because the terms Q F , Q ' F , and Q" contain curvature information 
about the surface of the moving body. For these reasons. Formulation 3 is used only 
for the transonic and supersonic regions of the moving body even though it is valid 
everywhere. 

COMPUTATIONAL STRATEGY 

In this section, the computational strategies for applying equations (2), (3), and 
(4) to propeller noise predictions are discussed. The section is composed of several 
subsections which highlight various analytical and computational facets of the solution 
procedure. 

Blade and Surface Pressure Modelling 

Calculations with ASSPIN begin by modelling user supplied blade geometry and 
loading data with high order, tensor product, least-squares splines. The modelling is 
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performed automatically by a modified version of the software appearing in reference 
12. Blade surface parameterization involves representing the Cartesian coordinates of 
blade points in terms of spanwise, ;/ 2 , and percentage chord, Q, curvilinear coordi- 
nates. Steady blade loads are parameterized in a similar fashion. The resulting two 
dimensional suiface functions arc easily evaluated and differentiated as required by the 
integrands of (2) and (3). 

Unsteady blade surface pressure due to nonaxial inflow is periodic with period 2n fu, 
where u: is the angular speed of the propeller. At each user defined blade point, unsteady 
pressure is input as a function of time over the entire period of revolution. To evaluate 
the unsteady pressure and its derivatives as required by ASSPIN, the Fourier time series 
of the unsteady blade surface pressure is computed. The Fourier coefficients, which are 
functions of the blade curvilinear coordinates ?/ 2 an d Q, are stored and interpolated as 
required. At present, ASSPIN computes the first five components of the Fourier series. 
If there are regions of time in which the surface pressure varies rapidly, for example if 
there is a pylon in front of the propeller, then additional Fourier components may be 
requited foi adequate representation. In these cases, ASSPIN can be easily modified 
to accommodate the increased resolution. 

Previous versions of the noise prediction code evaluated the unsteady surface 
pressure by straightforward three dimensional interpolation. Compared to steady 
surface pressure noise predictions, predictions in which unsteady loading was required 
used excessive amounts of computer time. This problem has been alleviated by the 
incorporation of the above strategy in which evaluation of the unsteady pressure has 
been reduced to a two dimensional interpolation problem. Unsteady surface pressure 
noise predictions now require only 1.25 times the execution time of steady surface 
pressure predictions. 

Coordinate Reference Frames 

There are three Cartesian reference frames used in the noise calculations. They 
are the medium fixed x-frame, the aircraft fixed X-frame, and the blade fixed 7 f-frame 
(Figure 3). Initially, it is assumed that the propeller center is at the origin of the 
x-frame. The x 3 axis lies in the direction of flight, with positive x 3 pointing in the 
direction of motion. The directions xq and x 2 are defined so that the x - frame is right- 
handed. The aircraft fixed frame is identical to the medium fixed frame initially and 
moves in the forward flight direction with speed V F . 

The origin of the blade fixed frame coincides with the origin of the aircraft fixed 
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frame at all times. The angle between the 77 3 axis and the X3 (or -Y3) axis is the inflow 
angle, a. The pitch change axis of the propeller blade always lies along the 772 ax i s > 
and r]i is chosen so that the if frame is right handed. Note that if the inflow angle is 
zero, then the 773 axis and the A'3 axis always coincide. For single rotating propellers 
and the forward propeller of CRPs, the blade is assumed to rotate counterclockwise in 
the 771772-plane with respect to an observer looking down the 773 axis. The aft propeller 
of CRPs is assumed to rotate clockwise. 



Figure 3. ASSPIN Reference Frames 


The mathematical relationships between the three reference frames are obtained 

by applying the appropriate translations and rotations. Let 77 denote the blade fixed 

— ¥ 

coordinates of an arbitrary blade point. At time r, let y(r) and Y(t) represent the 
locations of the moving blade point in medium fixed coordinates and aircraft fixed 
coordinates, respectively. If lo is the angular speed of the propeller, then the following 
relations exist: 

( 6 ) 



and 


Y(t) - A l A 2 {t)i% 


( 7 ) 
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where .4i and ^(r) are the 3x3 rotation matrices 


and 


cos a 0 — sin o 

Ai = I 0 1 0 

sin a 0 cos a 


COS LOT — sin cur 0 

A' 2 { T ) = ( sin cur coster 0 

0 0 1 


The velocity, V(r), of a blade point with respect to an observer in the medium fixed 
frame is obtained by differentiating (6) with respect to source time. Thus, 


V(r) 


dyij) 

dr 


= Ai 


d- 4 - 2 (r) 

dr 


V + 



( 8 ) 


where 


d^ 2 (r) 

dr 


— sin cur — costar 0 
cos cur —sin cur 0 

0 0 0 


In order for the integrands of (2) and (3) to be meaningful, all of the terms must 
be written in the same coordinate system. Most of the terms contain only blade 
geometric quantities. These quantities are easily computed in terms of blade fixed 
coordinates by evaluating the spline functions that were described in the previous 
subsection. Therefore, to facilitate the evaluation of the integrands, all terms that 
involve vector operations (e.g., dot products) are expressed in blade fixed coordinates. 
It is natural to write the Mach vector, A I, of a blade point and the moving observer 
point, ;r, in medium fixed coordinates. Thus, M and x must he expressed in blade fixed 
coordinates in order to maintain consistency. 

The Mach vector, expressed in terms of medium fixed coordinates, is obtained from 
(8). Since M is viewed as a free vector, it follows that its representation with respect to 
the aircraft fixed frame is the same as that of the medium fixed frame. Consequently, at 
time t, the Mach vector of a moving blade point, if, expressed in blade fixed coordinates 
is calculated by applying the inverse transformation of (7) to equation (8) which yields 


M { r )= W(r)A;' 9 -§A 



dA -2 (t) 

dr 


ij + A 2 VMi 1 



( 9 ) 
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where 


and 


a;'(t) = 



v = 


cos a 0 sin a 

0 1 0 

— sin cx 0 cos a 


Consider an observer point that moves with the aircraft. Suppose that the initial 
position of the observer with respect to the medium fixed frame is given by (xi , # 2 , x s)- 
At the observer time i, the observer point has translated V F t distance units in the £3 
direction. Therefore, the coordinates of the observer point at time t with respect to 
the aircraft fixed frame at the source time r are (ri, X 2 , £3 -f V F (t — r)). Thus, the 
location of the observer with respect to the blade fixed frame, ; r), is given by 


\(t; t ) = A 2 1 (t)A 1 



( 10 ) 


Let if be the blade fixed coordinates of a blade surface point, and let r be the source 
(emission or retarded) time, then, according to equation ( 10 ), the radiation vector, ? r , 
appearing in equations ( 2 ) and ( 3 ) can expressed in terms of blade fixed coordinates as 

r = x(t,T)-ff. (11) 

Retarded Time Calculation 

Computations with either Formulation 1 A or Formulation 3 involve the repeated 
calculation of retarded (or emission) times. For fixed observer position, r, and observer 
time, t, the retarded time, r, is a solution of the retarded time equation (RTE) 

|r| = \x(t) — y(r )| = c 0 (< — t) , ( 12 ) 

where y is a fixed source point whose emission time(s) are sought. Equation ( 12 ) states 
that an acoustic signal, traveling at the speed of sound, that leaves y at time r takes 
t — r time units to reach the observer position x at time t. It can be shown (reference 
13) that the RTE has exactly one solution for subsonic motion between source and 
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observer. However, if the relative motion is supersonic, then the RTE may have many 
solutions. At current design conditions, propfans operate at helical tip speeds of about 
1.15. It has been found that the RTE has at most three solutions at these speeds. 

For a typical supersonic helical tip speed case, retarded time calculations account 
for approximately 50% of ASSPIN’s total computational time. Furthermore, since the 
theoretical formulations require accurate retarded times to be useful, a robust algorithm 
for solving the RTE is necessary. In the ensuing discussion, the precise form of the 
RTE for propeller motion is developed, and an efficient method for finding its roots is 
presented. 

Suppose (?/i, r/ 2 , 7/3) are the blade fixed coordinates of a blade point and 
(a’ 1, X 2 , £3) are the initial coordinates, with respect to the medium fixed frame, of 
an observer point. If t is the observer time, then the RTE can be combined with 
equation (11) to yield 


r-r~ |xT<; r) - ff\ 2 = c 2 (t - r) 2 . 


(13) 


Let 


<f> = U! (t — t) , 


= x\ cos a + X3 sin a, 
x *2 = x 2 ? 

x 3 = * r 3 cos a ' — x l sin a, 

** = \Al 2 

V=y/ri f + vl 

x * 

i’x = tan -1 '-j, 
x \ 

and 

/ f -1 7 /2 

Wn = tan — , 

'/l 


then equation (13) can be written 


A(f> 2 + D(f> + C + cos ((f) + D) + E(f> sin (<f> + F) = 0, 


(14) 


where 


14 



c 2 - V 2 
A~° F 


B = 


2 rjx*uj' 2 ’ 

V F [.r* sin a + (*3 - 173) cos a] 


urjx 


x* 2 + ij 2 + (.T3 - 773) ] 


C = 

ior/x 

D =ip ,, - rp x + ut, 
Vr, sin a 


E = -^— 7 

U)X 

F =V’n + <X>t. 


and 


Equation (14) is an implicit, transcendental equation for 4> which cannot be 
solved analytically. Further complications arise when the motion of the source point 
toward the observer is supersonic. In this instance, (14) has multiple solutions and 
straightforward root finding algorithms, such as Newton’s method, must be augmented 
by analysis if all roots are to be resolved. A robust algorithm for finding all roots of 
the RTE has been developed by the authors and is discussed below. 


Since equation (14) was obtained by squaring equation (12), it follows that spurious 
solutions of (14) might exist. This problem is resolved by enforcing the causality 
condition, which states that a signal must be emitted before it can be observed. Thus, 
if 4>* is a solution of (14), then it must satisfy the relation 


sgn{u)(f)* < 0. 


(15) 


This is because cj> — u! (t — t), and T must be less than or equal to t if the causality 
condition is to be satisfied. 

Let g((f>) denote the left hand side of (14). Equation (15) gives limited information 
on where to search for roots of g. However, sharp bounds can be obtained by analyzing 
(14). It can be shown that 

cos ( (£> + D) + E<!) sin (<j> + F) = yEfy 2 + (1 + E 2 <j>f cos (<f> + D — t/>), 


where 


E\ — E cos 4'x , 


E -2 = E sin 4’x , and 




= tan 


Eil 

1 + e 2 4> 


15 



Therefore, (14) can be written as 


A<f) 2 + Bcj> + C 

~ r = = — COS ( <j) + D — V>). 

yjE 2 ^+(l + E 2 <j>) 2 


In absolute value, the right hand side of the above equation is bounded by one. Thus, 


Tb 2 + B 0 + C 

\J E 2 6 2 + (1 + E 2 4>) 2 


< 1 for all <j>. 


The roots of g that satisfy (15) are therefore bounded by the roots of 

{a.4> 2 + B<f) + C) = Ei4> 2 + (1 + E 2 (j>) 2 . (16) 

Equation (16) is a polynomial equation of degree four, whose roots are easily calculated. 
ASSPIN uses the computer program RPOLY, taken from reference 14, to solve (16). If 
rf'min anc ^ $ max are the two solutions of (15) and (16), and if 4>* n i Jl < <f >* nax , then any 
roots of g must lie in the interval <P* max \- 

The next step in the root finding process involves the approximation of g{<t>) 
° n timin’ ^ max ] by a polynomial of degree three. The function g is evaluated at 
four evenly spaced points on [<^* un , 4>max\ an< l fh e interpolating cubic polynomial 
calculated. Any extreme points on [<?*,„■„ , <^n tax ] of the cubic polynomial are then 
computed analytically. If extreme points exist, then they are used as initial guesses for 
finding the extreme points, if any, of g by Newton’s method. If extreme points of g are 
discovered, then the cubic polynomial is recalculated by replacing some of the original 
interpolation points by the extreme points. The advantage of this last step is discussed 
below. 

The final phase in the calculation of retarded times is to find the roots of the third 
degree polynomial analytically. Roots of the polynomial which lie in [<f>* nin , <f>* nax ] are 
then used as initial guesses to Newton’s method for determining the roots of g. 

It is difficult to distinguish between two roots that lie close together. Trying to 
resolve such roots using Newton’s method requires accurate guesses for the initial 
iterants. Otherwise, Newton’s method may diverge or keep converging to the same 
root. Two roots of a smooth function are separated by an extreme point. Therefore, 
using extreme points of g as interpolation points for the cubic polynomial enhances 
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the approximation of g in the neighborhood of the adjacent roots. Consequently, if the 
roots of the cubic near an extreme point of g are used as initial guesses, then Newton s 
method converges rapidly to the nearest roots of g. This technique ensuies that the 
search for all roots of the RTE is exhaustive. 

Spectral Analysis 

For a given observer point, x, translating with the propeller, the periodic acoustic 
signature generated by a propeller with N B blades is obtained by the superposition of 
the acoustic pressure signatures produced by the individual blades. If the blades are 
identical and symmetrically spaced and if p\{x, t) is the signature of one of the blades, 
then the cumulative pressure signal is given by 


A * / 

p'(x, t) = f + 

n= 1 '• 


2n(n — 1)\ 

y 


For each observer location, ASSPIN calculates p\ at a discrete number of equally spaced 
time points, and applies the above formula to obtain the combined acoustic signature. 

Before proceeding to the details of how is calculated, the spectral anal}/ sis of p is 
considered. Acoustic output from ASSPIN consists of discrete pressure signatures for 
the loading, thickness, and combined noise, and the sound pressure levels as a function 
of harmonic number for each noise component. The temporal Fourier series of the 
acoustic pressure can be written 


p {x , t) = ^ a n (F) cos uj n t + b n (x)s\nu n ti 
n — 0 

where ui n = nN B io. The Fourier coefficients are given by 

T 

a n (x) = | t p'(x, t) cos wntdt, (17a) 

0 

T 

b n (x) = J Pi*' t)smu> n tdt, (17b) 

0 


and 


where T 


2?r 

TV? 
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Since p* is known for a discrete number of observer time points, equations (17a) 
and (17b) must be evaluated numerically. ASSPIN computes the Fourier coefficients 
by Simpson’s rule (reference 11). Simpson’s rule provides a higher degree of accuracy 
than Fast Fourier Transform (FFT) methods and is not limited to a specific number 
of time points ( viz. for some k ). As a result, greater flexibility and a higher degree 

of accuracy is attained relative to FFT methods. It is also worth mentioning here that 
the additional computational time required by Simpson’s rule is negligible relative to 
the entire noise prediction process. 


The n-th sound pressure level for the total acoustic pressure, SPL n , is defined by 
the equation 

SPL n = 10 log, 0 - " (x) + h ^\ 

^re f 

where p re j = 2 x 10 5 Pascals is the reference pressure. Sound pressure levels for the 
thickness and loading components of the noise are defined similarly. 


Discretization and Numerical Integration 

The blade suiface is separated into two regions - the subsonic root region and the 
transonic/supersonic tip region. The spanwise location 


* 

12 


sgu(u>) 


UJ 


,2 


— cos 2 

r 


a — V F sin 


“)’ 


where c* = (1 - e)c 0 with e a small user supplied parameter and a is the inflow angle, 
serves as the boundary separating the two regions. This separation point insures that 
the high speed portion of the blade, which might require the use of Formulation 3 to 
evaluate the noise contribution, is adequately disjoined from the subsonic region. Both 
regions are then subdivided into equally spaced spanwise and percent chord intervals 
according to user specifications which results in a panelled discretization (Figure 4a). 

For a specified observer point, x, the one blade periodic acoustic pressure signal, 
P\{ x i t), is calculated at a discrete number of observer time points comprising the 
period of revolution for one blade. For this purpose, an adaptive observer time grid is 
employed. The adaptive grid is a function of observer location, operating conditions, 
and source geometry^. To facilitate the spectral analysis of the acoustic signals, the 
results for the adaptive observer time grid arc interpolated onto an evenly spaced 
observer time grid at the end of the calculations. 

The idea behind the adaptive grid strategy is to determine the range of observer 
times (if such a range exists) in which blade points have multiple emission times. For 


18 



(a) Blade Discretization 
N1SUB N1SUP 


intervals intervals 



Figure 4. ASSPIN Computational Grids 

observer times within this range, the panel contains transonic and/or supersonic sources 
requiring the use of Formulation 3. Panel sources are subsonic for observer times outside 
of the range, and the acoustic pressure for the panel at these observer times can be 
computed with Formulation 1A. 

The range of observer times in which a fixed observer/source point pair has multiple 
emission times can be found by analyzing equation (14). Again denote the left hand 
side of (14) by the function g(<f>). The coefficients D and F of (14) are functions of 
observer time, therefore the solution(s), <j), of g(<f>) = 0 is a function of t. A typical 
example of this functional dependence is shown below in Figure 5. In this illustration, 
a situation is depicted in which the RTE has multiple solutions when t 6 [t m i n , t mar } ■ 
From Figure 5, observe that 

-* oo at t = t m i n and t = t max . (IS) 

at 

Differentiating (14) with respect to t yields 

d(j> sin(</> + -D) — E<fr cos(<ft + F) 

dt ~ U 2 A<f> + D - sm(<f> + D) + E sin(<£ + F) + E<f> cos (<f> + F) 
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Figure 5. Multiple Emission Time Illustration 


Therefore, it follows from (18) that 


2A.<f> -f- B — sin (0 + Z?) + E sin (0 T F ) + E<f> cos (0 T F) = 0 (19) 


at t — tmin tmax- 

The observer times t mjn and tmax are obtained by solving equations (14) and (19) 
simultaneously for t and 0. It can be shown that the observer time can be eliminated 
from (14) and (19) resulting in the equation 


[w + sin0v) ^.40 2 + B(j) + C^j — cosiJ> x ( 2A0 + B) j 
+ ^(cos ifrx ~ E) (yA(f) 2 + B<f> + C’j + (E(f> + sin-0 x ) (2A<f> + B) 


COS 0a; (COS 0j; — E) + (£0 + sin 0a;) 2 


( 20 ) 


Equation (20) is a sixth degree polynomial in 0 and ASSPIN finds its roots with 
RPOLY. The solutions of (20) that do not satisfy (15) are discarded and the observer 
times t m i n and t max are then found with either (14) or (19). This procedure is repeated 
for many source points on the blade. The final adaptive observer time grid for the blade 
is based on the minimum t mtn and maximum t max of all source points tested. 

The acoustic pressure, when viewed as a function of observer time, has steep 
gradients in regions of observer time in which the blade sources emit signals that 
approach the observer at transonic or supersonic speed. For this reason, the period 
of time in which the blade contains points that have multiple emission times is finely 
discretized. Conversely, the observer time period in which the blade exhibits subsonic 
behavior is coarsely discretized. 

The adaptive observer time grid strategy is a novel feature for high speed propeller 
noise prediction in the time domain. Its implementation in ASSPIN has reduced 


20 



computer execution time by approximately 50 percent relative to previous versions 
of the noise prediction codes. 

Consider a fixed panel on the discretized propeller blade. For a fixed observer point 
x, the acoustic signal generated by the source panel is computed by looping through the 
adaptive observer time grid. For each observer time point of the adaptive grid, ASSPIN 
determines if the noise contribution for the panel is to be computed by Formulation 
1A or Formulation 3. If the condition 


(1 — M r ) ret < e, e prescribed 

is satisfied at some source point, then Formulation 3 is used for the panel at that 
observer time. Switching back and forth between the formulations is performed 
automatically. This procedure is repeated for each blade panel and the resulting 
acoustic signals summed. 

If Formulation 1A is to be used, then the integrals of (2a) and (2b) are computed 
via Gauss-Legendre (G-L) quadrature (reference 11). A user prescribed number of 
G-L nodes and weights are calculated as part of ASSPIN preprocessing. The nodes 
are distributed over the panel (Figure 4b), and the integrands of (2a) and (2b) are 
evaluated at the nodes and summed according to G-L integration requirements. 

Due to the complex nature of the integrands of (3a) and (3b), a simpler numerical 
quadrature scheme is employed for calculating Formulation 3 integrals. The F curves 
of equation (4) and the acoustic planform boundaries are approximated by piecewise 
straight lines and the integrals of (3a) and (3b) calculated by the trapezoidal rule 
(reference 11). Increased resolution is obtained by a further subdivision of the panel 
(see Figure 4c) according to user specifications. 

The evaluation of the Formulation 3 surface integrals by the collapsing sphere 
method is a multiphase process. First, for a fixed observer time, the retarded times 
are computed at each vertex of the spatially refined panel. For each subdivided panel 
(refer to Figure 4c) the retarded times at the vertices are sorted in ascending order. 
Thus, the period of time in which the collapsing sphere intersects a particular panel 
is established. This period of source time is then discretized and the F curve, z.e., 
the intersection of the panel and the collapsing sphere, determined for each discretized 
source time point. The position of the F curve within a panel is calculated by testing 
the panel edges for intersection with the collapsing sphere. If an edge tests positive, 
then the intersection point is approximated by linear interpolation and after a second 
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edge is found, the location of the F curve is complete. The two intersection points are 
connected by a line segment and the integrands evaluated at the endpoints. Geometric 
quantities, such as normal vectors and curvatures, are linearly interpolated based on 
their values at the panel vertices. 

Counter-Rotating Predictions 

There are several restrictions with regard to counter-rotating propeller (CRP) noise 
predictions that must be emphasized. It is assumed that both the forward and aft 
propellers operate at the same angular speed and opposite in direction, have the same 
number of blades, and their pitch change axes (PCA) are aligned initially (Figure 6). 



Figure 6. Counter-Rotating Propeller Configuration 


These restrictions are imposed to facilitate the spectral analysis of CRP acoustic 
pressure waveforms. If nonlinear acoustic interactions between the propellers are 
neglected, then the resultant CRP signal is simply the linear superposition of the 
individual signals of the two propellers. It is assumed that aerodynamic interactions 
between the two propellers are accounted for via user supplied blade surface pressure. If 
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the user is interested in CR.P predictions where these assumptions cannot be satisfied, 
then the signals from the individual propellers can be obtained by operating ASSPIN 
in the single-rotation mode and adding them together in the desired fashion external 
to ASSPIN. In some instances, ASSPIN can be easily modified to accommodate special 
circumstances . 


A schematic representation of the noise prediction process is shown in Figure 7. 



Figure 7. ASSPIN Flowchart 


PROGRAM DESCRIPTION 

Details of program usage, input specifications, and program output are described 
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in this section. 


Operating Instructions 

User input consists of a one line identifier, two groups of NAMELIST variables 
which define operating conditions and grid sizes, blade geometry data, loading data, 
and observer coordinates. The structure and content of the input file will be discussed 
in the ensuing subsections. The user input file resides on FORTRAN logical unit 
number 5. 

Standard program output is written to FORTRAN logical unit number 6. Addi- 
tional output files containing intermediate results must be opened by the user in the 
calling program. These files are associated with FORTRAN logical unit numbers 7 
and S. ASSPIN output consists of a banner page, a one line identifier, NAMELIST 
parameters, computational grid information, computed power and thrust coefficients, 
thickness, loading, and combined noise spectra and signals. 

To access ASSPIN, the user’s calling program must have the FORTRAN statement 

CALL ASSPIN( X, NXDIM ). 

The variable X is a one dimensional dynamic storage array whose dimension, NXDIM, 
must be at least 

max [3 * NOBS + 9 * NT + 20 * N1S * N2S + 3 * inax(NlS,N2S) + 7 * (IGAUSS) 2 + 1, 
NEPTS * NQPTS * (5 + 2 * NTPTS) + 6 * NEPTS + 2 * NQPTS + NTPTS], 

The capitalized quantities appearing in the above expression are ASSPIN input 
variables and are defined in the next subsection. 

Input Specifications 

User input is divided into four categories: 1) physical characteristics of the propeller 
configuration and ambient medium, 2) grid size parameters for the various numerical 
schemes, 3) blade geometxy and either steady or periodically unsteady blade surface 
pressure data, and 4) observer coordinates. Details of all user supplied parameters and 
data are provided below. 

NAMELIST PHYSCAL Parameters. 

The physical characteristics of the propeller and undisturbed medium are input to 

ASSPIN via the FORTRAN statement NAMELIST PHYSCAL. In Table 1 below, 
the NAMELIST PHYSCAL parameters are described. The FORTRAN variable names 
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appear on tlie left hand side of the table and the variable’s definition, type, and units 
(if applicable) are on the right side of the table. 

Table 1. Input Parameter Description - NAMELIST PHYSCAL 


VARIABLE 

DESCRIPTION 

BTIP 

Real. Distance, along the PC A, from the propeller axis to the 
propeller tip in meters. 

BINNER 

Real. Distance, along the PCA, from the propeller axis to the 
propeller root section in meters. 

VF 

Real. Forward flight velocity of propeller in meters per second. 

RPM 

Real. Propeller angular velocity in revolutions per minute. 

CO 

Real. Sound speed of undisturbed fluid in meters per second. 

RHOO 

Real. Density of undisturbed medium in kilograms per cubic meter. 

BETA 75 

Real. Blade angle at 3/4 span in degrees. 

NB 

Integer. Number of blades. 

NOBS 

Integer. Number of observers. 

STEADY 

Logical. STEADY = TRUE => Steady surface pressure. 
STEADY = FALSE => Unsteady surface pressure. 

XINFLOW 

Real. Inflow angle in degrees (sec Figure 3). 

NCRP 

Integer. NCRP - 0 => SRP prediction. 

NCRP = 1 => Front propeller of CRP prediction. 
NCRP = 2 =t> Aft propeller of CRP prediction. 

PROPDIS 

Real. Distance between CRP propellers in meters. 


NAMELIST GRID Parameters. 

The computational grid sizes are defined by the user through the FORTRAN 
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statement NAMELIST GRID. Grid parameters are defined and their default values and 
restrictions, if any, are given in Table 2. The effect of varying the grid size parameters 
has been studied and the results reported in reference 15. 

Table 2. Input Parameter Description - NAMELIST GRID 


VARIABLE 

DESCRIPTION 

NT 

Integer. Number of observer time points for single blade pressure 
signature. Restriction: NT > 2. Default: NT = 125 * NB. 

NS 

Integer. Number of sound pressure levels to be output. 
Default: NS = 30. 

N1SUB 

Integer. Number of equally spaced spanwise intervals for the subsonic 
portion of the blade (Figure 4a). Restriction: N1SUB > 1. 

Default: N1SUB = 2. 

NlSUP 

Integer. Number of equally spaced spanwise intervals for the 
supersonic portion of the blade (Figure 4a). Restriction: NlSUP > 1. 
Default: NlSUP = 3. 

N2SUB 

Integer. Number of equally spaced percentage chord intervals for the 
subsonic portion of the blade (Figure 4a). Restriction: N2SUB > 1. 
Default: N2SUB = 4. 

N2SUP 

Integer. Number of equally spaced percentage chord intervals for the 
supersonic portion of the blade (Figure 4a). Restriction: N2SUP > 1. 
Default: N2SUP = 4. 

NlS 

Integer. Number of equally spaced spanwise points for Formulation 3 
panels (Figure 4c). Restriction: NlS > 2. Default: NlS = 10. 

N2S 

Integer. Number of equally spaced percentage chord points for 
Formulation 3 panels (Figure 4c). Restriction: N2S > 2. 
Default: N2S = 10. 

EPSILON 

Real. If EPSILON > 1 — M r for some point on a panel then 
Formulation 3 is used. Default: EPSILON = 0.05 

IGAUSS 

Integer. Number of Gauss- Legendre nodes used to compute 
Formulation 1A integrals (Figure 4b). Restriction: IGAUSS > 1. 
Default: IGAUSS = 7. 
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Blade Geometry and Surface Pressure. 

Blade geometry is input by specifying the shape and orientation, with respect to 
the blade fixed reference frame, of airfoil sections. The number of sections to be input, 
NEPTS, is determined by the user. 

Airfoil sections are generated by cutting the PC A with planes perpendicular to the 
PCA. The orientation of a section is determined by specifying its distance along the 
PCA from the propeller axis, 172, the local chord length, CH, the blade fixed coordinates 
of the leading edge, rjf E and 7/3 £ , and the difference, A/ 3 , between the local blade 
angle, /?, and the blade angle at 75 % span, / 3 75 . The above measurements are assumed 
dimensional with distance measured in meters and angles measured in degrees. Note 
that the innermost cutting plane must be at the spanwise location t ] 2 = BINNER and 
the outermost cutting plane is at i] 2 — BTIP. These quantities are shown schematically 
in Figure 8. 


r\2 fixed 


T )2 fixed 



(a) Configuration for SRP and forward 
propeller of CRP 


(b) Configuration for aft 
propeller of CRP 


Figure 8. Airfoil Orientation Specifications 

The upper airfoil surface is defined by specifying the distance from the chord line 
to the upper airfoil surface, y i; , as a function of relative distance along the chord, 
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Q. By convention, Q — 0 corresponds to the leading edge of the section and Q — 1 the 
trailing edge. Also by convention, points lying above the chord line are assigned positive 
distances, while those lying below the chord line are assigned negative distances. The 
lower airfoil surface is defined in the same manner. Distance between the chord line 
and the lower surface is denoted by y L (see Figure 9). The quantities y v and y L are 
assumed to be nondimensionalized by the user with respect to local chord length. They 
are input at a discrete number of relative chord points, {Qi}^' PTS , which are the same 
for each section. The leading edge point and the trailing edge point must be included. 



Figure 9. Airfoil Coordinate Specifications 


Blade surface pressure must be supplied by the user at each ( 772 , Q) point for which 
blade data is input. The user must nondimensionalize surface pressure by the quantity 
p 0 u> 2 (BTIP) 2 . For each spanwise cut, steady surface pressure on both the upper and 
lower blade surfaces, p u and p L respectively, are input as functions of Q. 

Unsteady surface pressure is assumed to be periodic with period 2n. The user 
is responsible for this temporal nondimensionalization. For each blade point, the 
surface pressure is input as a function of time. Discretized time points, {U}^.T PTS , 
comprising the period are input by the user and must be identical for each blade point. 
Nondimensionalization of the unsteady surface pressure is the same as for the steady 
case. 

Descriptions of the geometric and surface pressure input parameters are summarized 
in Table 3. 
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Table 3. Geometry and Loading Input Parameter Description 


VARIABLE 

DESCRIPTION 

NEPTS 

Integer. Number of spanwise data points (airfoil sections). 
Recommend: 15-30. 

NQPTS 

Integer. Number of relative chord data points. 
Recommend: 15-50. 

NTPTS 

Integer. Number of time points for unsteady pressure data. 
If STEADY = TRUE, then NTPTS - 1. 
if STEADY = FALSE, then NTPTS > 12. 

Recommend: 15-50. 

Q 

Real. Relative chord data points. 
Nondimensionalizcd by CH. 

Q 1 = 0 and Qkqpts = 

m 

Real. Spanwise data points in meters. 

DINNER < 772 < BTIP. 

CH 

Real. Local chord length in meters. 

A/? 

Real. Difference, in degrees, between local blade angle and /? 75 . 

vf E 

Real. rj\ coordinate of leading edge in meters. 

n 3 LE 

Real. 773 coordinate of leading edge in meters. 

Vl 

Real. Distance between chord line and lower airfoil surface. 
Nondimensionalizcd by CH. 

Vu 

Real. Distance between chord line and upper airfoil surface. 
Nondimensionalizcd by CH. 

Vl 

Real. Pressure on lower airfoil surface. 
Nondimensionalizcd by /? 0 u? 2 (BTIP) 2 . 

Pu 

Real. Pressure on upper airfoil surface. 
Nondimensionalizcd by p 0 ^ 2 (BTIP) 2 . 

t 

Real. Time points for unsteady surface pressure. 
Nondimensionalizcd by fj = 0 and t NTPTS = 2n. 
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Observer Coordinates. 

The acoustic observers are assumed to translate with the aircraft. Their coordinates 
are expressed in meters and are specified by the user with respect to the aircraft fixed 
reference frame. 

Input File Structure 

The order and format of the user supplied input file is discussed in this subsection. 
Both SRP and CRP predictions are considered. Examples involving both cases appear 
in the next section. 

SRP Predictions. 

1. 80 character identifier - FORMAT(1X,A80). 

2. NAMELIST PHYSCAL - NCRP = 0 for SRP Predictions. 

3. NAMELIST GRID. 

4. Blade coordinates and surface pressure data. 

4a. NEPTS, NQPTS, NTPTS - F0RMAT(1X, 3110). 

4b. Relative chord grid - F0RMAT(1X, E15.7). 

Qi = 0.0 

QnQPTS = 1-0 

4c. t? 2 , CH , A/3, t)£ e , ri£ E - F0RMAT(1X, 5E15.7). 
y L (Qi), y v {Qi), p L (Qi), Pu(Qi) - format(ix, 4E15.7). 


y l(Q NQPTS )> Vu (Q NQPTS )> Pl(Q NQPTS )> Pu^Q NQPTS ) ' F0RMAT(1X, 4E15.7). 

The quantities p L and p v refer to steady blade pressure. If the blade loads are 
unsteady, then ASSPIN does not use p L and p Lr . However, properly formatted 
numbers must appear in these locations. 

Item 4c. is repeated for each of the NEPTS airfoil sections. 

4d. Time point grid for unsteady surface pressure - F0RMAT(1X, 5E13.6). 

READ(5)(T(K), Iv = 1, NTPTS) 
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T(l) = 0.0, T(NTPTS) = 2tt. 

If STEADY = TRUE, then 4cl. is omitted. 

4e. Unsteady surface pressure - F0RMAT(1X, 5E13.6). 

DO K = 1,NTPTS 
DO I = 1,NEPTS 

READ(5)(p l (I,J,K), J = 1,NQPTS) 

READ(5)(p a (I,J,K), J = 1,NQPTS) 

If STEADY = TRUE, then 4e. is omitted. 

5. Observer coordinates - Unformatted. 

READ(5,*)((f(I,J), J = 1,3), I = l,NOBS) 

af(I,J) is the J-th Cartesian component, with respect to the aircraft fixed 
frame, of the I- til observer. 

CRP Predictions. 

1. 80 character identifier - FORMAT(ASO). 

2. NAMELIST PHYSCAL for forward propeller - NCRP = 1. 

3. NAMELIST GRID for forward propeller. 

4. Blade coordinates and surface pressure data for forward propeller. 

See format and instructions for SRP predictions. 

5. Observer coordinates. 

See format and instructions for SRP predictions. 

The origin of the aircraft fixed frame is assumed to be at the hub center 
of the forward propeller (see Figure 6). 

6. NAMELIST PHYSCAL for aft propeller - NCRP = 2. 

NOBS and NB are the same for both propellers. RPM for the aft propeller must be 
the negative of RPM for the forward propeller. 

7. NAMELIST GRID for aft propeller. 

NT must be the same for both propellers. 
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8. Blade coordinates and surface pressure data for aft propeller. 

See format and instructions for SRP predictions. 

9. Observer coordinates. 

See format and instructions for SRP predictions. 

This information must be the same as for the forward propeller. 

EXAMPLES 

Three sample calculations that demonstrate the capabilities of ASSPIN are pre- 
sented in this section. In each example, the operating conditions and acoustic re- 
sults are displayed. The source code and input and output files for the three cases 
can be obtained upon request. The authors can be reached at the email address 
mhd@sabre.larc.nasa.gov. 

SRP Prediction with Axial Inflow 

In this example, the noise generated by the SR-7L propeller (reference 16) was 
computed. The propeller is eight bladed, has a diameter of nine feet, and is designed 
to cruise at Mach 0.8. The NAMELIST PHYSCAL are as follows: 

BTIP = 1.3716, 

BINNER = 0.3429, 

VF = 239.4, 

RPM = 1679.0, 

CO = 299.3, 

RHOO = 0.372, 

BETA75 = 57.0, 

NB = 8, 

NOBS = 1, 

STEADY = TRUE, 

XINFLOW = 0.0, 

NCRP = 0, 

PROPDIS = 0.0. 

The acoustic observer is located 0.62 diameters from the propeller tip and 0.25 
diameters aft of the propeller plane. The coordinates of the observer with respect to 
the aircraft fixed frame are (0.0, 3.06, -0.686). The NAMELIST GRID parameters are 
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set to their default values (see Table 1). Aerodynamic loads were computed with the 
Adamczyk code (reference 17). 

The computed thrust and power coefficients for this example are C T = 0.53 and 
C p = 1.85, respectively. Acoustic waveforms and spectra for the thickness, loading, 
and combined noise are shown in Figure 10. 

SRP Prediction with Nonaxial Inflow 

The SR-7L propeller was again used for this case. The propeller was set at an inflow 
angle of 4.6 degrees and the unsteady blade loads were computed with the Whitfield 
code (reference 18). The NAMELIST PHYSCAL parameters are as follows: 

BTIP = 1.3716, 

B INNER = 0.3429, 

VF = 239.4, 

RPM = 1679.0, 

CO = 299.3, 

RHOO = 0.372, 

BETA75 = 57.8, 

NB - 8, 

NOBS = 1, 

STEADY = FALSE, 

XINFLOW = 4.6, 

NCRP = 0, 

PROPDIS = 0.0. 

The acoustic observer is the same as in the previous case. ASSPIN grid parameters 
were again set to their default values. 

The computed thrust and power coefficients for this case are C T = 0.43 and 
C p = 1.75, respectively. Acoustic waveforms and spectra for the noise components 
are displayed in Figure 11. 

CRP Prediction with Axial Inflow 

The results computed here were produced by a model counter-rotation propeller 
system. The blade loading, which includes the aerodynamic interaction, was obtained 
from Adamczyk’s code. For the forward propeller, the NAMELIST PHYSCAL param- 
eters are as follows: 
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BTIP = 0.3064, 

DINNER = 0.137, 

VF = 242.9, 

RPM = 8283.0, 

CO = 338.3, 

RHOO = 1.039, 

BETA 75 = 56.3, 

NB = 8, 

NOBS = 1, 

STEADY = TRUE, 

XINFLOW = 0.0, 

NCRP = 1, 

PROPDIS = 0.1037, 

and for the aft propeller: 

BTIP = 0.2973, 

BINNER = 0.1322, 

VF = 242.9, 

RPM = -8283.0, 

CO = 338.3, 

RHOO = 1.039, 

BETA 75 = 53.9, 

NB = 8, 

NOBS = 1, 

STEADY = TRUE, 

XINFLOW = 0.0, 

NCRP = 2, 

PROPDIS = 0.0. 

The acoustic observer is located 1.62 diameters from the forward propeller tip and 
slightly aft of the propeller plane of the forward propeller. The coordinates of the 
observer with respect to the aircraft fixed frame are (0.8, 0.0, -0.05). 

The computed thrust and power coefficients for the forward propeller are C T = 0.34 
and C p = 1.25, and for the aft propeller C T = 0.39 and C p = 1.36. Acoustic results for 
this case appear in Figure 12. 
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Figure 10. Acoustic Results for SRP with Axial Inflow 
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Figure 11. Acoustic Results for SRP with Nonaxial Inflow 
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Figure 12. Acoustic Results for CRP with Axial Inflow 
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CONCLUSIONS 


The computer code ASSPIN predicts the noise produced by advanced technology 
propellers in either the single- rotation or counter-rotation mode. It is based on two 
theoretical formulations of Farassat, which are valid in the near and far fields and for 
axial and nonaxial inflow conditions. One formulation is valid for subsonic helical tip 
speeds and involves integrals over the surface of a propeller blade. The other is used for 
transonic and supersonic helical tip speeds and is written as surface and line integrals 
over the acoustic planform, generated by the blade, and its boundary. 

Efficient, computational algorithms are used in ASSPIN to evaluate the theoretical 
formulations numerically. Blade geometry and surface pressure data are approximated 
by high order splines. Accurate numerical integration schemes and an adaptive 
observer time grid have been developed to reduce computer execution time. The 
most computationally intensive portion of ASSPIN involves the accurate calculation 
of retarded times. In a typical noise prediction, the retarded time equation must be 
solved about 500,000 times and in some instances has multiple solutions. ASSPIN 
employs an exhaustive root finding procedure which accurately calculates the retarded 
times in a few iterations. 

Details on the use of ASSPIN are given in this paper. Input to ASSPIN consists 
of operating conditions, observer coordinates, and blade geometry and surface pressure 
data. ASSPIN produces acoustic signatures and spectra of thickness, loading, and 
overall noise. 
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